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Abstract 

Sea ice breaks up and regenerates rapidly during winter conditions in the Arctic. 
Analyzing satellite data from the Kara Sea, we find that the average ice floe size 
depends on weather conditions. Nevertheless, the frequency of floes of size A is a 
power law, N ~ where r = 1.6 it 0.2, for A less than approximately 100 

km?. This scale- invariant behaviour suggests a competition between fracture due 
to strains in the ice field and refreezing of the fractures. A cellular model for this 
process gives results consistent with observations. 
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Polar ice is a dynamic phenomenon linked to sea currents and global climate. 
Even in the coldest winter conditions, the ice fields in the Arctic seas are un- 
stable to fragmentation and drift. Cracks in the polar ice have been observed 
spanning the entire polar region. Here we study the fragmentation process 
within a small section of the Arctic ice in order to obtain an empirical de- 
scription of fracture and healing of sea ice. 
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Fig. 1. Left: One of the regions in the Kara Sea where the sateUite data have 
been collected, delimited by the square. Right: Identification of ice floes from two 
successive pictures of the above area, which were taken in the period 24-27 February 
1994.(c)ESA. For details, see text. 

Many fragmentation processes are governed by the irreversible formation of 
cracks[l,2,3], typically modeled by simplified approaches to the dynamical in- 
terplay between stress redistribution and crack formation[l,4,5,6,9] across the 
fragmenting system. Ice fields in the Arctic seas represent a different situation 
as emerging cracks in addition may undergo rapid healing due to refreezing. 
Hence, in these systems, a steady-state situation ensue where the fragmen- 
tation and the healing process balance each other. Note, however, that even 
though the two processes are opposite, they do not lead to the system being 
reversible with respect to time: Fragmentation and healing are very different 
processes. 

In winter conditions, Arctic sea ice typically freezes to a thickness of approx- 
imately two meters. Ice fields are subject to stresses due to wind and sea 
currents. Such stresses cause cracks in the ice field, and further fragmenta- 
tion occurs when ice floes collide with each other. Ice floes merge when the 
exposed water between floes freezes, and the cracks heal. The competition 
between fracture and healing results in a distribution of ice floe sizes that 
depends on weather conditions. 

In the right panel of Fig. 1 we show a 10000 km^ ice fleld in the region of the 
Kara sea depicted on the left side of the picture. The ice floes were identifled 
from synthetic aperture radar (SAR) pictures taken with 3-day intervals in 
the period January-March 1994 by the European Space Agency (ESA) ERS-1 
satellite. The ice fioes were determined by measuring the displacements of a 
regular grid of points 500 m apart over a 3-day interval. Thin arrows on the 
right panel of Fig. 1 show displacements during the period 24-27 February 
1994. From this set of up to 40,000 ice displacements rigid areas are identified 
by searching for vector regions which transform through rigid transformations 
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Fig. 2. Frequency of ice floes at one particular time (24 February, 1994), squares, 
and for the time averaged frequency of ice floes as a function of their area A, circles, 
calculated from satellite data. Also included are calculations with the simple model 
of ice fracture and refreezing described in the text, full line. 

{i.e. translations followed by rotations) [10,11]. Each of these rigid areas is 
identified as an ice floe. 

The data allow us to determine 29 different ice field patterns from time series 
of satellite data over the entire period. The ice field is quantified by counting 
the number of ice floes of area A occurring in each sample. The distribution 
of ice fioes, illustrated in Fig. 2 for the case of the same data set as in Fig. 1, 
shows scaling with ~ A^'^ and r = 1.6 ±0.2. A similar behaviour was found 
for all 29 sets, with variations in the width of the scaling regime from frame 
to frame. 

The temporal variation in ice field patterns can be quantified by, for exam- 
ple, the variation in average ice floe size. This variation reflects the widely 
different weather conditions, with wind speeds reaching up to 30 m/s, and air 
temperatures ranging between -10 and -40 C. 

The overall distribution, represented by the circles in fig. 2, obtained by sum- 
ming the whole data set, shows a power-law behaviour that extends over 
more than two orders of magnitude, with a similar exponent as the individual 
frames. 

Over the last decade, there has been a strong interest in systems with compet- 
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Fig. 3. Rules for cracking and healing depending on the local geometry of randomly 
selected pre-existing cracks. The circle represents the selected point. See text for 
details. 



ing mechanisms and irreversibility. The present system is of this type. Such 
systems may organize themselves to a steady state characterized by scale in- 
variance, as is the case for turbulence [7] or systems exhibiting self organized 
criticality [8]. 

The appearance of scaling behaviour in the ice floe distribution may be de- 
scribed by a simple model where fracture and healing compete to shape the 
ice field. A feature of ice fioe dynamics is that stresses build up between 
fragments, making both crack healing and crack formation initiate along pre- 
existing cracks. Our dynamic model for fracture and healing is defined on a 
two-dimensional grid. Depending on geometry, we either heal the crack or cre- 
ate a new one perpendicular to it, by applying the following rules at each time 
step: 

i - A random point on a boundary between two fragments is selected, which 
defines the interacting region at the current time step. 

ii - As shown in fig. 3, a new crack always appears when: 

(iia) the overlapping region is less than the smallest touching fragment's 
side; 

(iib) the chosen point is situated on the corner of an ice fioe. 
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In both cases, one of the two cracking possibihties is selected with equal 
probability. 

iii - Healing happens in either of the cases displayed on fig. 3, and stated 
below: 

(iiia) the interacting region entirely contains the touching border of one of 

the floes; 

(iiib) the selected point corresponds to the intersection of three ice floes. 
In this case, one of the three healing possibilities illustrated is randomly 
chosen. 

iv - Either cracking or healing may occur with equal probability if the selected 
point is found in particular configurations, such as: 

(iva) the interacting region of one of the fragments is "L" shaped; 
(ivb) both fragments are of the same lateral size and are in full lateral 
contact. 

As displayed on fig. 3, different cracking configurations are randomly chosen. 
V - Nothing is done if the point corresponds to the place where 4 cracks 
meet. 

When a new crack is created, it starts at the selected point and propagates 
until an other crack is found, or the border of the mesh is reached. A crack 
may only start at the mesh borders if there are no cracks in the system. This 
happens in the first steps of the evolution since we start with a single ice floe 
which spans the whole mesh. If a crack is healed, the line corresponding to the 
crack is followed and erased until no open line is left, such as in the healing 
case shown on panel (iVa) of fig. 3. 

As a result we obtain a dynamical steady-state behaviour of the ice floe field, 
where ice fioes of all sizes appear. The distribTition of fragment sizes is a power 
law with r ~ 1.5, as is depicted by the full line in fig. 2. When a large mesh 
is employed, the system exhibits a huge number of possible configurations. 
Although, as already mentioned, we start the dynamics with a very particular 
configuration, irreversibility is observed after a large number of steps because 
tracing back the original configuration would be virtually impossible follow- 
ing the rules of the model. Therefore, after the transient regime, the system 
attains a dynamically changing state which bears no fingerprint of the initial 
conditions, or of any other instant configuration in its remote past. 

We remark that no explicit attempt has been made to relate crack formation 
to neighboring fioe sizes with a preference for fragmenting the largest, or the 
thinnest of them. Neither does the model take into account the possible pref- 
erence for equal size fragments to induce cracks, as suggested by [4]. We stress 
that when a crack is induced as illustrated in fig. 3, one of the floes is chosen 
with equal probability. When refreezing occurs the crack to be healed is se- 
lected in the same way. The fragmentation process, however, assumes that all 
cracks initiate on pre-existing one dimensional fractures, a feature reminiscent 
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of the boundary induced fragmentation discussed by [6]. Our geometric rules 
intend only to incorporate some basic physical properties in the model. For 
instance, rule (ii-a) is an attempt to take into account the torsion along the 
touching fragments due to the center of mass displacement of the ice floes, 
which would very likely induce new cracks. Energy sharing among the frag- 
ments would tend to make new fracture development less probable when three 
ice floes meet. Instead, the splash of cold water is more likely to occur, leading 
to refreezing of the cracks as implemented in rule (iii-b). Similar qualitatively 
considerations apply to the other geometries. Instead of introducing ad-hoc 
weights to cracking or healing in each case, we simply chose the most proba- 
ble case. We have checked that a characteristic size scale appears if the rules 
are inverted, whereas only finite size effects seem to be observed when one 
adopts those presented above. 

The present model is related to a variety of models that have been employed to 
describe crack formation in different systems [1,5,6]. It differs by considering 
an interaction proportional to contact between fragments, and, in particular, 
by proposing a mechanism that naturally balance crack formation and healing. 

In summa, the ice floe system presents us with a case study for a largely 
unexplored class of fracturing/healing phenomena with steady state dynamics 
— phenomena that would include continental plate dynamics on time scales 
of hundreds millions of years. 
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